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The tensile behavior of a thin web-like paper network was simulated by two simple mathematical models. The 
mesh distortion, drop in tensile force and energy loss resulting from breakage of a network junction were 
calculated. These results were used to formulate two parameters for characterizing interfiber adhesion: a parameter 
averaging the network energy losses incurred in a series of bond breaks when the network is elongated and a 
parameter averaging the force drops. The effect of mesh size, local bond adhesive force, and size and shape of the 
specimen network were calculated. These results based on model studies were used to interpret behavior observed 
in an actual paper network. 
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1 . Introduction 

There is a need to develop parameters that can be used to characterize quantitatively the adhesion between 
fibers in a sheet of paper. Such parameters would be useful in selecting, blending and treating batches of 
pulp to produce a paper with improved interfiber adhesion. This could facilitate the recycling of paper pulps, 
leading to important savings in energy and materials. 

As a result of some recent research [1—3] a new technique is evolving for obtaining these parameters. A 
handsheet in the form of a very low-density web is made from the pulp to be evaluated. If a specimen from 
this handsheet is elongated to break in a sensitive tensile tester, a force-elongation curve containing numerous 
jags is obtained. It is assumed that each jag is caused by the breakage of a bond between fibers constituting 
the handsheet network. Thus if the values of the change in force corresponding to each jag are averaged, a 
force drop parameter characteristic of interfiber adhesion might be obtained. The energy lost by the network 
as the result of a bond break can be found by integrating under the force-elongation curve in the vicinity of a 
jag. An average of these energy losses could also be used as a measure of bond adhesion between two fibers. 

The selection of these parameters was an intuitive process. Their applicability must be judged on the basis 
of experimental evidence. Moreover, what the parameters actually measure is not well understood. Some 
questions that arise are: How is the strain energy distributed throughout the network after a bond break? How 
are the energy loss and force drop related to the local force at the bond just before break? What is the effect 
of mesh size? What is the optimum specimen size and shape for tensile tests to determine the characteristic 
energy? 

In the following, two crude mathematical models of a fibrous network are described, and their behavior 
under elongation calculated. The concepts developed and the results obtained add to an intuitive 
understanding of the force drop and energy parameters. 



2. The Parallel-Spring Model 

The parallel-spring model is depicted schematically in figure 1. It consists of a number of springs in a 
series-parallel arrangement. The dots indicate where the springs are bonded to each other and to a series of 



* This work was sponsored by the U. S. Department of Knergy. 

1 Figures in brackets indicate literature references at the end of this paper. 



299 



I 

C\J 



B 



2M-I 



FIGURE 1. Schematic representation of parallel- spring model. 

Model has '2M — 1 parallel columns of springs attached to 2/V — 1 rigid transverse bars. 



rigid bars. There are 2M — 1 junction points on each bar, and there are 2/V — 1 bars. In the unstrained state 
each spring has a length /, called the mesh length. The spring constant is given by k/l. The model is 
elongated so that each spring is extended an amount /8, where 5 is the strain. The bond B then breaks, and 
the two springs previously joined at that point are inactivated. 
The initial force sustained by the model before break is 



The initial elongation of the model is 



F b = (2M - l)k8 



M = 2(/V - 1)18 



The initial energy stored in the model is 

E b = V 2 (2M - I)k8-2(N - 1)18 
= 2(N - 1){2M - l)l- 1 /2kd 2 



(1.1) 



(1.2) 



(1.3) 



After the break the elongation of the model is maintained constant, and the force drops to a value F a . The 
horizontal bars to which the springs are attached remain horizontal. In the two segments of the model on 
either side of the broken bond the force F a is sustained by 2M — 2 parallel springs, so that 



F a = (2M - 2)k 



(AL), 
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(AL) ' = (2F^ (L4) 



where (AL) t is the elongation of the segment. 

The elongation in the other 2(/V-2) segments is given by 



FJ 
(AL) 2 = , * ^, (1.5) 



(2M - l)/c 



where (AL) 2 is the elongation of the segment. 
The total elongation of all the segments is 



2(/V -2)FJ 2FJ 

2(7V - 1)/8 = -7 —?- + 



Thu 



(2A# - 1)A: (2M - 2)A 



, /V - 2 2M - 1 

(2M - 1)A:8 = F b = 4- — -\F a 

1 TV — 1 2(M- 1 (/V- 1)| " 



2fl/ - 1 1_ 

l+ 2(M-l)(/V- 1) TV - if" 



2(M- 1)(/V - 1) 



He 



F fl = TF b (1.6) 



vhere 



2(M-1)(A^-1) 

T = — 1.7) 

2(M - l)(N - 1) 4- 1 

T is the fraction of the force remaining after a bond break. The force drop resulting from the break is given by 

(1.8) 













AF = 


F b 


~T) 




The 


energy in the 


two 


segments 


of the model adjoining 1 


he bn 


►ken bond is 












1 /2(A£) 1 -F fl 


= Vi 


Fa 

(2M - 


H 
-2)k 


and th< 


» energy 


in th( 


• <>t 


her 2(/V - 


2) segments is 


















l h(M) 2 -F a 


= Vi 


Fa 2 l 

(2M - l)k 














301 







Thus the total energy is 



F 2 1 F 2 I 

E a = 2(N - 2)-V2 - a - + 2- V2 - a - 
; (2M - 1)* (2M - 2)* 

/V - 2 1 1 F 2 l 



2 + -V2 

l2M- 1 2M- 21 A; 



2iV- 1) [yV- 2 2M - 1 / - / 

= - -\ +- : — T 2 - (2M- l) 2 k 2 8 2 

2M - 1 I N - 1 2(M - 1)(7V - 1) I 2k 

= 2(N - l)(2M- l)/-i-f 2 -V»*6 2 



£ a = 2(7V - 1)(2M - l)/f- V2 A:6 2 

= TE b (1.9) 

In order to describe the way strain energy is distributed in a model it is useful to associate energy with each 
of the junction points in such a way that the sum of the associated energies is equal to the total energy stored. 
At any given junction point the associated energy is equal to one-half the sum of the energies in the springs 
joined at that point. The factor of one-half is introduced because when the associated energies are summed 
over all the junction points, the energy of each spring is counted twice. At junction points on the top and 
bottom boundaries there is only one spring, but the factor of one-half is still necessary. 

A dimensionless average associated energy E can be defined by means of the relation, 

E a = (2M - 1)(27V - 1)E- 1 l2kb 2 l (1.10) 

The factor (2M — 1)(2/V — 1) gives the total number of points in the network. The factor 1 /2k8 2 l is the 
energy stored in one spring before a bond break, and is introduced into the expression to make E 
dimensionless. By equating relations (1.9) and (1.10) an alternative relation, 

2N - 2 - 

E = T (1.11 

2N - 1 

is obtained. 

The energy loss resulting from a bond break is, from eq (1.9), 

l±E = E b - E a = (1 - f)E b (1.12) 

An alternative form is obtained from eqs (1.3) and (1.10). 

27V — 1 - 
AE = 2(N - 1)(2M- 1)/(1 £)-V2A;§ 2 (1.13) 

3. The Square-Network Model 

The square-network model is depicted schematically in figure 2. It consists of springs forming a square 
network. There are 2M — 1 columns of junction points and 2N — 1 rows of junction points. The depicted 
network has 7 columns (M = 4) and 7 rows (N = 4). The unstrained length of each spring is /, and the spring 
constant is given by k/l. The model is elongated by an amount 2(/V — l)/8. The bond at point B then breaks, 
and the four springs previously joined at that point are inactivated. The resulting configuration is calculated 
by a computer program described below. 
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FIGURE 2. Schematic representation of square- net work model. 

Model has C 2M — 1 parallel columns of springs. Each column has 2/V - 1 junction points which, excepting the end points, are 
attached to transverse springs. 



In order to discuss the model it is convenient to identify junction points by the coordinates i = 0, 1, "•• m; 
j = 0, 1, ••• n. The corresponding numbers/ = 1, 2, ••• M; J = 1, 2, ••• N, however, are retained for use in 
the formulas for force and energy, as they are more convenient for counting purposes. Because of symmetry it 
is only necessary to calculate the configuration of the upper right-hand quadrant; i.e., the junction points 
corresponding to positive values oft andy. 

When the model is strained each of the junction points moves from its original position by an amount £,j / 
in the horizontal direction and 77^ / in the vertical direction. Thus when the model is first extended and no 
bond break has occurred, 






(2.1) 



After the bond break the values of fy, Tj/j readjust except for certain constraints at the boundary and the 
i = 0,y = axes of symmetry. These are 



£;„ = 0. 

Vin = id 

Vm = 0. 

it.-i = in 
Vi,-i = Vn 

toi = 0. 



(2.2) 
(2.3) 
(2.4) 
(2.5) 
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b-U blj 

V-u = Vv 



(2.6) 



The values of £, r\ attained after bond break are found by solving a system of equations in which the 
horizontal and vertical forces at each junction point are balanced out, plus the boundary conditions (2.2-2.6). 
To formulate the equations applicable at the point i, j , let 



r a = Vo + A^) 2 + a^ 2 



R B = VA B f + (1 + A B 7,) 2 



R c = VA r £ 2 + (1 + A c t,) 2 



(2.7) 



R D = V(l +A^) 2 + A,rf / 

where A^ is a difference operator involving the indices (i, i — 1), A B the indices (j, j — 1), A c the indices 
(j + 1,/) and A/) the indices (i + 1, i). Thus, for instance, A A £ means £ tj — £i-ij. The quantities/?/ are the 
distances from a point ij to an adjacent point. The horizontal force balance equation is then 

k(R D - 1)(1 + A^ + k(R c - 1)A^ = k(R B - 1)A B £ + k(R A - 1)(1 + A^) 



R D 



Rc 



Rb 



(1 + A D £)(1 - \/R D ) + A r £(l -l/Rc) - A B £(1 - \/R B ) - (1 + A fl £)(l - \/R D ) = 0. (2.8) 
Similarly for the vertical components 

A d t}(1 - \/R D ) + (1 + A c tj)(1 - l/« f ) - (1 + A b t,)(1 - \/R B ) - A^(l - \/R A ) = 0. (2.9) 

Equations (2.8) and (2.9) apply at the general junction point ij 9 but when i = 0, m orj = 0, n they must be 
modified to conform with boundary or symmetry conditions. 

It should be noted that the horizontal forces acting at a point i, j are dependent mostly upon the values of 
?i-i,i» £u an d fi+i,j« Therefore eq (2.8) which has the form 

f(€i-ij,m-i,ii f«. Vij, fi+u. T?i+i,j, €u-i>yu-u &,j+u Vij+i) = °- ( 2 - 10 ) 

was expressed as a Taylor's series expansion in terms of only these three most important variables to become 

where the Afy means fy ( * +1) ~ fy (fc) and A^-j j and A^ l+1 j have similar meanings. / (fc) , 1 etc. are 

evaluated as functions of i; (k \ r\ {k) . The values of £ (fc+1) are obtained by solving the system of eqs (2.11) 
corresponding to different values of j . 

A set of five of these equations, for instance, would look like 



BjX ] + C,X 2 


= *i 


A 2 X 1 + B 2 X2 + C2X3 


= ^2 


A zx 2 + B 3 x 3 + C3X4 


= K 3 


A 4X3+ B 4 x 4 + C 4 % 5 


= K 4 


/4 5 * 5 + 5 5* 5 


= K S 
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This set of equations is easily solved by a Gauss-Jordan reduction. 

Equations (2.11) are solved to obtain values of £ {k+l \ Similarly eq (2.9) can be expanded in terms of 
Tfjj^-i , Tjij and Tfaj+i • The set of equations obtained for different values of i can be solved to obtain values of 



„<*+« - i d f\ i (1 f 



^(/f+n j ne va j ues () { J I — — 1 — — j n e q (2.11) for £ and the corresponding equation for 7] are then 

\MijJ \dyu/ 

evaluated using g ik+] \ ^k+i) an( j t ^ e Drocec l ure iterated to obtain g (k+2 \ r) ik+2) . 

During the calculation process solutions for £ and r) may be obtained such that R A , R B , R c , or R D becomes 
less than 1 in value. This means that there is a spring in compression. This situation is not permitted as it is 
more realistic to assume that the spring would buckle rather than resist a compressive force. Therefore if there 
is an R less than 1, modified equations in which the influence of the spring is removed are used in an 
alternate calculation. 

The configuration of the network can be calculated in other ways. For instance eqs (2.8) and (2.9) could be 
expanded in the variables fy, r\ y. Improved approximations for £, r\ could be obtained by solving only two 
simultaneous equations at each point, but a very large number of iterations would be required to obtain very 
precise values of £, r\. On the other hand the equations could be expanded in terms of all 10 variables. In this 
way a system of approximately MN equations in MN unknowns would be obtained. If these equations had a 
numerically stable solution and round-off errors were not too severe, only a few iterations might be required. 

The method of solution that was used was a compromise between these two extremes. It appears to be 
stable but in most cases requires 1(X) or 2(H) iterations to obtain satisfactorily precise values. When situations 
are calculated in which the initial strain 8 is of the order of one percent, more than 200 iterations may be 
required. In this situation the values of some Rs are very close to unity and very small values of R — 1 are 
obtained. Thus it might be better in this case to rewrite eqs (2.8) and (2.9), expressing terms of the form 
(1 — l/R) as a truncated series, although this was not done for the calculations reported here. 

After sufficiently accurate values of £ and r\ were obtained, the tensile force T in at the top of each column 
of springs was calculated using the formula 



T m = k(\ +nS -i?(.»-,)0 - I/Rb) 

R B = Vff.,,-1 + (1 + "8 - Vi,n-lf 
When i = 0, the value of T 0n simplifies to 



(2.12) 



The force after break F n then becomes 



T 0n = k(n8 - T) .„-i) (2-13) 



F a = T 0n + 2 2 T in (2. 14) 



The fraction of the force remaining after a bond break T is given by the ratio F a /F b . 

The average associated energy E is found from its definition, eq (1.10), after calculating the associated 
energy for each junction point and summing over all the junction points to obtain the stored energy after break 
E a . 

The value of E can be obtained by a different calculation, if F a is known as a function of 5 from a series of 
calculations on a given model. The energy after break E a is found by integrating under the curve of F a versus 
initial elongation A L = 2(/V - l)/o\ Thus, 

^2(A/-l)/fi 

£ a = F a -d[2(N- 1)18] 

.'o 

(2.15) 

-8 



= 2{N - 1)(2M -\)kl \ tddd 
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From the definition eq (1. 10) E is then given by 

2/V 



E = 



2 2 



2/V - 1 8 2 Jo 



T8dd 



(2.16) 



In the square-network model T is not a constant independent of 8 as it is in the parallel-spring model. If T 
were constant, eq (2. 16) would reduce to the relation eq (1.11) given for the parallel-spring model. 

It should be noted that the behavior of the square-network model is the same as that of the parallel-spring 
model before the central bond breaks. Therefore eqs (1.1), (1.2) and (1.3) are the same for both models. The 
quantity T is defined by eq (1.6) and the quantity E by eq (1.10) in both models. Therefore eqs (1.8) and 
(1.13) are valid for both models. Relation (2.16) is a general definition valid for both models, but relations 
(1.7) and (1.11) apply to the parallel-spring model only. Therefore eqs (1.9) and (1.12) also are valid only for 
the parallel-spring model. 



4. Square-Network Model Calculations 

A number of calculations were carried out for different square-network models to determine the number of 
iterations required for accurate values of T and E. The model for which M = II, N = 21, I = 1 was studied 
the most extensively. Results of calculations for this model are given in table 1. These results indicate that for 

Table 1. Values of T and E as Functions of Number of Iterations for a Square-Network 
Model M = 11, N = 21,1 = 1. 



8 


No. of 
Iterations 


f 


E 


0.01 


100 


0.95227 


0.92923 




200 


.95254 


.92923 


.03 


100 


.95779 


.93112 




200 


.95714 


.93100 


.05 


20 


.96738 


.93606 




50 


.96679 


.93583 




100 


.96614 


.93567 




200 


.96564 


.93561 


.10 


20 


.97806 


.94643 




50 


.97829 


.94638 




100 


.97827 


.94637 


.15 


20 


.98316 


.95270 




50 


.98330 


.95263 




100 


.98332 


.95262 


.20 


10 


.93787 


.97979 




20 


.98371 


.95679 




50 


.98594 


.95644 




100 


.98595 


.95641 




200 


.98597 


.95640 



Table 2. Values of T and Efor Some Square-Network Models 8 = 0.20, 1 = 1, 
200 Iterations. 



M 


N 


f 


E 


11 


11 


0.97554 


0.92306 


11 


21 


.98597 


.95640 


21 


21 


.99311 


.96639 


11 


41 


.99238 


.97634 


21 


41 


.99626 


.98221 
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values of 8 greater than 0.05, LOO or 200 iterations are required to obtain values of T and E accurate to more 
than four significant digits. When 8 is 0.05 or less, more than 200 iterations are desirable. Check calculations 
were also carried out on the other models studied, and 200 iterations were determined as sufficient. 

Calculations were carried out on five different models. The initial strain 8 for these models was 0.20. A 
mesh length of/ = 1 was assumed (or the calculations. Values of M and N were varied in order to study how 
T and E depended on the shape and on the number of junctions in the model. Results of the calculations for 
these models are given in table 2. 

The configuration of the model for which M = 11, TV = 21 is plotted in figure 3. Only the upper right-hand 
quadrant of the model is shown, as the complete model is symmetrical with respect to this quadrant. The 
distortion due to the bond break is seen to be concentrated along the central column of springs (/ = 1). 
Buckling occurs (R < 1) in the regions where the springs are indicated by dotted lines. Although not easily 
discernable in the figure, the central width (/ = 1) has been reduced to 0.956 times the original width, as a 
result of bond break. 
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FIGURE 3. Configuration of the positive quadrant 
of a square-network model for which M = 11, 
N — 21, after a central bond break. 

The broken bond is designated by the open circle in the lower left 
corner. Springs in which buckling occurs are designated by dashed 
lines. 
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Figure 4. Energy map of positive quadrant of a square-network model for which M — 11, N — 21, after a central bond break. 

Vertical-axis scale gives ratio of associated energy at a junction point after bond break, divided by the associated energy before break. 



Table 3. Values of T and E as Functions of Initial Strain for a Square-Network Model 
M = 11, N =21, I = 1. 



8 


t 


E 


E 

eq(2.16) 


0.01 


0.95254 


0.92923 


0.92931 


.03 


.95714 


.93100 


.93230 


.05 


.96564 


.93561 


.93658 


.08 


.97477 


.94270 


.94328 


.10 


.97827 


.94637 


.94683 


.13 


.98173 


.95052 


.95078 


.15 


.98332 


.95262 


.95273 


.18 


.98507 


.95510 


.95503 


.20 


.98597 


.95640 


.95626 



Note: T for the parallel-spring model is 0.99751 

E for the parallel-spring model is 0.97318 
Note: Values for °= 0.01, 0.03, 0.05 and 0.20 were obtained after 200 iterations. All 

other values were obtained after 100 iterations. 



An energy map of the model for which M = 11, AT = 21 is shown as an isometric plot in figure 4. As before 
only the upper right-hand quadrant of the model is shown. The dimensionless value plotted at each junction 
point is the ratio of the associated energy after bond break to the associated energy before bond break. 
According to this map the associated energy remains constant at most of the junction points, it increases at 
points immediately to the right and left of the broken bond and decreases to a lower level at junction points 
above and below the broken bond. 
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FIGURE 5. Configuration of the positive quadrant 
of a square- network model for which M -21, 
N = 41, after a central bond break. 

The broken bond is designated by the open circle in the lower left 
corner. Springs in which buckling occurs are designated by dashed 
lines. 



Figure 5 shows the configuration after bond break, and figure 6 shows the energy map of the model for 
which M = 21, N = 41. This model, of the same shape as the previous model, has approximately four times 
as many junction points. The corresponding figures for the two models are similar. 

Table 3 presents values of t and E as functions of initial strain 8 in the case of the model for which 
M = I], /V = 21. A force-elongation curve for this model after one central bond has been broken can be 
plotted using data from table 3. The force has the value (2M - 1)7*8, and elongationthe value 2(N " l)/8. 
The essential characteristics of the curve are shown by figure 7, which is a plot of f 8 vs 8 for this model. 
Unfortunately the values of t are so close to unity that this curve lies very close to the force-elongation curve, 
shown as a dashed line, of the model before the central bond was broken, and the nonlinear character of the 
curve is not readily apparent. 

The non-linear character of the force-elongation curve can be demonstrated by plotting the difference in 
force before and after break as a function of 8. This is shown in figure 8, in which (1 - f)d is plotted vs 8. As 
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Figure 6. Energy map of positive quadrant of a square-network model for which M = 21, N = 41, after a central bond break. 

Vertical-axis scale gives ratio of associated energy at a junction point after bond break, divided by the associated energy before break. 
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FIGURE 7. Force- elongation curve after central bond break in a square- network model for which M = 11, N =21. 

T8 is plotted vs 8, where 8 is the elongation per unit length. Dashed line shows force-elongation curve before bond break. 
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Figure 8. Curve of force drop due to bond break vs elongation of a square-network model for which Ml — II , /V 

(1 — T)8 is plotted vs 8, where 8 is the elongation per unit length. Dashed line shows force drop for equivalent parallel-spring model. 
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a matter of interest (1 — T) for the parallel-spring model (M = 11, /V = 21) was calculated. Its value was 
0.0024938. The linear funetion 0.00249388 is shown as a dashed line on the plot. The value of T for the 
parallel-spring model exeeeds the values of T calculated for the square-network model. Evidently the foree 
drop incurred by a bond break in the square-network model is greater than the force drop incurred in the 
parallel-spring model. 

The area under the T8 curve in figure 7 is equal to / T8d8, and according to eq (2.16) the average 



associated E is equal to 



2/V-2 2 



times this quantity. Thus values of E can be calculated from T and 8 



27V - 1 8 2 

according to eq (2. 16) and the result compared with the values of E in table 3, which were calculated from the 
associated energies present after break. This should provide a check on the accuracy of the computer 
calculations. The results of this calculation are given in the last column of table 3. The values of E obtained 
by the two methods are seen to be in fair agreement. 



5. Discussion of Results 

7. Configurational Distortions and the Redistribution of Assoeiated Energies as the Result of a Bond Break. 

It has been assumed that in a low-density paper network of the type tested experimentally [1-3], in most 
cases a bond break results only in local distortions, and most of the energy loss is also concentrated locally. 
It would be desirable to demonstrate this with a mathematical model, but the square-network model is not 
suitable for this purpose. The distortions and energy losses are concentrated along the column of springs 
directly in line with the broken bond. Possibly a better model would be one in which the meshes were of a 
different shape, such as a hexagon, and so arranged that load bearing fibers could not be easily aligned in the 
direction of stretch. Such a model would require less stretching force, have a nonlinear force-elongation 
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curve, and probably have a greater tendency to buckle laterally. It would better simulate the behavior of a 
thin paper network, but at the price of greater mathematical complexity. 

In the energy maps, figures 4 and 6, there is an increase in associated energy at the junctions to either side 
of the broken bond. These bonds therefore should break immediately after breakage of the central bond, and 
bonds farther out should break successively in a tearing action that proceeds across the model network. This 
tearing action is probably common to many models, and actually occurs in real paper networks. In the paper 
networks however it seldom goes to completion, because the density of fibers is not sufficiently uniform. 

2. The Effect of Mesh Length. 

In experimental work specimens of standard dimensions (usually 2X1 cm) are elongated, and the applied 
force recorded as a function of the elongation. Whenever a bond breaks the force decreases sharply, and the 
energy loss resulting from the bond break can be determined by integrating under the force-elongation curve 
in the vicinity of the force drop. Certain of these energy losses are averaged to obtain a parameter 
characterizing bond strength. 

If it is assumed that all of the breaks occur at approximately the same local force level, and the specimens 
have uniform mesh size, this energy parameter might provide a reliable indication of the relative bond 
strength. For instance bonding could be studied in thin handsheets of the same mass per unit area made from 
a pulp that had been subjected to various beating treatments. However if one attempts to compare bonding in 
standard handsheets made from different pulps, a difficulty arises. One pulp may be coarser than the other. 
Thus a handsheet of a standard mass per unit area made from a coarse pulp would have fewer bonded 
junctions and larger meshes on the average than would a handsheet made from a finer pulp. In order to 
compare energy parameters for these different pulps it is necessary to know what is the effect of mesh size. 

To study this effect in models, consider two parallel-spring models: model A; M = 11, N = 21, / = 1 and 

model B;M = 21, N = 41, / = 0.5. These models have the same external dimensions, but the mesh length 

of model B is half that of model A. Let the two models be extended the same amount, so that 8 and therefore 

the local force k8 is the same in each. From eqs (1.1), (1.7) and (1.8) the force drop in model A is 21/c8/401 

~ 0.05A;8, and in model B is 4U8/1601 ~ 0.026&8. From eqs (1.3), (1.7) and (1.12) the energy loss in 

840 1640 

model A is Va k8 2 ~ 2. 1-Va *8 f , and in model B is Va k8 2 ~ 1.02'Va k8 2 . Evidently model B 

401 1601 

with half the mesh size of model A also has approximately half the force drop and half the energy loss. 

In any parallel-spring model, from eqs (1.1), (1.7) and (1.8) 

AF= w - m . 2(¥-D + i 

2(M - l)(N - 1) + 1 2(M - l)(N - 1) + 1 



which for large M becomes 



AF ~ - 1 — k8 

N - 1 



Substituting the initial length L = 2(A^ — 1)/ into this expression gives 



21 
AF~-k8 (4.1) 



From eqs (1.3), (1.7) and (1.12) 



2(7V - 1)(2M - 1)/ 1 . 
AE = — ~ r~ — k8 2 

2(M - 1)(7V - 1) 4- 1 2 

2(N - l)(2M - 2) + 2(N - 1) 1 
_ 1 . _ ^ 

2(M - 1)(7V - 1) + 1 2 
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which for large M becomes 



AE ~2l'V2k8 2 



(4.2) 



According to eqs (4. 1) and (4.2), if any two parallel-spring models for which M is large are extended so that 
the local force k 8 is the same in each model, the ratio of the energy losses AFJAF2 in the two models will he 
the same as the ratios of the respective mesh lengths /1//2. If in addition the initial lengths of the two models 
are the same, i.e., 2(7Vj — 1)/! = 2(/V 2 — 1)^ tne ratios of the force drops will ho equal to /Jl 2 - 

Table 4 gives the results of calculations for a number of parallel-spring models, all of which have the same 
initial length. According to these results the force drops and energy losses are almost linearly proportional to 
the mesh length /, despite the relatively low value of M in some of these models. It is interesting to observe 
that values of total force before break F b and the total energy before break F b do not influence these results 
significantly. The values of F b and E b are proportional to 2M — 1 however, so this is just a manifestation of 
the near independence of the values of AF and AE on the values of A/, for M sufficiently large 
(e.g.,M= 11). 

Table 5 gives results of calculations for some square-network models having the same initial length. The 
M, /V and / parameters for these models are the same as those in table 4. The force drops and energy losses for 
these models decrease as the mesh length / decreases, but are not linearly proportional to the mesh length. 
For instance the three models for which TV/ = 11 have / ratios of 1.0:0.50:0.25, AF ratios of 1.0:0.57:0.31 
and AF ratios of 1.0:0.64:0.37. The two models for which M = 21 have/ ratios of 1.0:0.50, AF ratios of 1.0: 
0.54 and AF ratios of 1.0:0.58. Apparently the dependence of AF and AF upon / becomes more linear as M 
increases. 

Results for the parallel-spring and square-network models suggest that the force drops and energy losses 
observed in a paper fiber network have an approximately linear dependence on a characteristic mesh length. 
The actual dependence, however, must be determined experimentally. 

The effect of mesh length can also be deduced by a different argument as follows: Select a model (/If, N) for 
which calculations have been made so that T and E are known for a given x; i.e., T = T(8), F = F(8). Select 
a mesh length /. The initial length of the model then is 2(/V — 1)/. For each mesh length selected the model is 
elongated an amount AL = 2(/V — l)/6\ This assures that the local tension at the bond break is always the 
same(k8). From eqs (1.1), (1.8) and (1.13) 



AF = (2M - 1)(1 - t)kb 



AE = 2(/V- 1)(2M- 1)/ 1 - 



2N 



2N 



2/2 



(4.3) 
(4.4) 



TABLE 4. Force Drops AF and Energy Losses AE for Some Parallel-Spring Models. 



Initial length^ = 2(/V - 1)/ = 20. for all models. 





M 


N 


/ 


Fjkb 


E b /0.3k8 2 


AF/k8 


AE/0.5k& 




11 


11 


1.0 


21. 


420. 


0.1045 


2.090 




I' 


21 


0.5 


21. 


420. 


.0524 


1.047 




21 


21 


0.5 


41. 


820. 


.0512 


1.024 




1 1 


41 


.25 


21. 


420. 


.0262 


0.524 




21 


41 


.25 


41. 


820. 


.0256 


.512 





Table . c 


. Force Drops AF and Energy Losses 


AE for Some Square- 


Network Models. a 




M 


N 


I 


F b /k8 


EJO.Sktf 


AF/k8 


AE/0.5k& 


11 


11 


1.0 


21. 


420. 


0.5137 


12.931 


11 


21 


0.5 


21. 


420. 


.2946 


8.270 


21 


21 


0.5 


41. 


820. 


.2825 


7.749 


1 1 


41 


0.25 


21. 


420. 


.1600 


4.811 


21 


41 


0.25 


21. 


820. 


.1533 


4.520 



Initial length L = 2(/V - 1)/ = 20. for all models. 
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As T and E are functions of 8 which is kept constant, AF is a constant independent of mesh length and AF 
is directly proportional to the mesh length. This conclusion is valid for either the parallel-spring or square- 
network type of model. 

The conclusion just obtained seems to conflict with that found by the analysis presented above. However in 
that analysis a number of models of different mesh length but having the same initial length were elongated, 
whereas in the present analysis the initial length of the model depends on the mesh length. 

This latter deduction may be useful in experimental work, as illustrated in the following situation: Suppose 
that two handsheets have been manufactured from the same pulp stock. Let the first handsheet have twice the 
areal density as the second, so that its characteristic mesh length, is one-half as large. According to the 
analysis above, bond breaks in both handsheets should occur at the same force level, but a specimen from the 
second handsheet must be extended twice as much as one from the first in order to attain this force. Thus the 
average energy loss per bond break in the second handsheet should be twice that in the first. This suggests 
that energy parameters found by tests on handsheets of different areal density can be scaled to find the value 
corresponding to a standard mesh length. 

3. The Effect of Sample Size and Shape 

If several models, each with the same mesh length /, are elongated to the same tension, the force drops and 
energy losses are given by eqs (4.3) and (4.4). For square-network models it is necessary to know the 
appropriate values of T and E before the effects of sample size can be predicted. For parallel-spring models 
however, simple expressions for T and E are known, so that for large values of M eqs (4.1) and (4.2) apply. 
Equation (4.1) predicts that the product of the force drop by the length of the specimen, LAF, is a constant 
for all models having the same mesh length /, and eq (4.2) then states that the energy losses are constant. 

Table 6 gives values of the force drop and energy loss for some parallel-spring models. Note that the 
quantities LAF/k8l and AE/0.5k8 2 l have the same value in this table because of the choice of units. It is 
apparent that AF and AE for parallel-spring models are closely predicted by eqs (4.1) and (4.2). The three 
models for which M = 11 (and width W = 2(M — 1)1 = 20/) have values for LAF and AE of approximately 
2.095, and the two models for which M = 21 (W = 40/) have values of approximately 2.048. Evidently LAF 
and AE approach a value of 2.0 with increasing M, as they should according to eqs (4.1) and (4.2). 

Table 7 gives values of the force drop and energy loss for some square-network models. The values of LAF 
and AE for these models are roughly constant, but there seems to be some dependence upon both M and N. 
For instance the three models for which M = 11 (W = 20/) have lengths that increase in the ratio 20:40:80. 
The LAF values for these models are 10.27:11.79:12.80. The two models for which M = 21 (W = 40/) have 
length ratios of 40:80 and LAF values of 11.30:12.27. Similar results are obtained for the AE values. 





Table 6. Values 


of Force Drop AF and Energy Loss AE for Some Parallel-Spring Models. 


M 


N 


L/l 


AF/k8 


LAF /k 81 


AE/0.5k8 2 l 


11 


11 


20 


0. 1045 


2.090 


2.090 


11 


21 


40 


.0524 


2.095 


2.095 


21 


21 


40 


.0512 


2.047 


2.047 


11 


41 


80 


.0262 


2.097 


2.097 


21 


41 


80 


.0256 


2.049 


2.049 



Table 7. Values of Force Drop AF and Energy Loss AE for Some Square-Network Models. 



M 


n 


L/l 


AF/k8 


LAF /k 81 


AE/0.5k8H 


11 


n 


20 


0.5137 


10.27 


12.93 


11 


21 


40 


.2946 


11.79 


16.54 


21 


21 


40 


.2825 


11.30 


15.50 


11 


41 


80 


.1600 


12.80 


19.25 


21 


41 


80 


.1533 


12.27 


18.08 
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Table 8. Values of &E/AF for Some Square-Network Models. 



M 


N 


HI 


1 A£ 
18 AF 


2 AF 
L8AF 


11 


11 


20 


12.59 


1.259 


11 


21 


40 


28.07 


L.403 


21 


21 


40 


27.43 


1.372 


1 1 


41 


80 


60.14 


1.503 


21 


41 


80 


58.96 


1.474 



The results for the square-network models suggest that various values of LAF and AF obtained by tests on 
speeimens of the same characteristic mesh length but of different size and shape, are roughly comparable. For 
most accurate results however all test specimens should have' the same standardized dimensions. 

For the parallel-spring models, from eqs (1.1), (1.3), (1.8) and (1.12) 



A£_£&_ 2(/V - l)/8 L8 
AF ~ F b ~ 2 2 



(4.5) 



This equation states that for the parallel-spring models the ratio AF/AF is proportional to the average 
elongation at which bond breaks occur, regardless of the value of M. Experimentally it is desirable to have 
significantly large force drops for a given energy loss, so the ratio AF/AF should be kept small by making/, 
as small as feasible. 

For the square-network models, from eqs (4.3) and (4.4) 



AF 
AF 



2/V 



I 



27V - 2 



F)/(\ - T) 



L8 
2 



(4.6) 



Table 8 gives values for the ratio AF/AF for some square-network models. Parallel-spring models have 
values of AF/AF that are proportional to the initial length L in accordance with eq (4.S), but for the square- 
network models the dependence of AF/AF upon L is not quite linear, and there also seems to be a small 

2 [AE\ 
dependence upon M. This is best shown by comparing values of the quantity — - — 1 , which for parallel- 



L8 \ AF) 

spring models is equal to unity but for the square-network models in table 8 varies between 1.2 and 1.5 
depending upon the model. 

Of most importance, however, is the confirmation that the quantity AE/AF can be kept low by choosing a 
small value of L. This situation is limited somewhat in experimental tests. A specimen length at least twice 
the specimen width is preferred, in order to avoid excessive stress distortion near the clamps. 

Although the results of tables 7 and 8 indicate that the values of LAF, AE and AF/LAF are insensitive to 
the value of M, there are circumstances where it is advantageous to vary the width of the specimen, as in the 
following experimental situation: Suppose that tests are being conducted on a given specimen, but the force 
drops observed are small, so that they cannot be measured with much certainty. The calculation of energy loss 
requires that the force drop be known as accurately as possible. Thus the experimental length should be small 
in order that the force drop be large. Further improvement is then achieved by decreasing the specimen 
width. This decreases the value of force at break F b without significantly affecting the value of the force drop. 
The amplification of the recording instrument can then be increased and the force drop better resolved. 



4. Dependence of Energy Loss Upon the Local Breaking Force. 



27V 



1 



According to eq (4.4) the energy loss AF is proportional to the product of the quantity I 1 

by the square of the local breaking force k8. For parallel-spring models E and therefore the quantity 

/ 2/V - 1 -\ . . . 

I 1 — E I is independent of the value of this force. Thus for these models energy loss is proportional 
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LOCAL FORCE , ARBITRARY UNITS 

FIGURE 9. Curve of energy loss due to bond break vs breaking force for a square-network model 

for which M = 11, N =21. 
2N - 1 - . 

( 1 E )8 is plotted vs elongation per unit length S, which is proportional to bond breaking force. Dashed line is plot 

27V - 2 
of K8 2 vs 8, and dotted line a plot of A '8' 5 vs 8. The constants K and K ' are adjusted so that ordinates of the plots are equal at 
8 = 0.10. 



to the square of the local force at break. For square-network models however the value of E depends upon the 

value of 8, so this simple quadratic relationship no longer applies. 

/ 27V — 1 -\ , 
Figure 9 is a plot of the quantity I 1 — — E 8 2 as a function of 8. The ordinate in this case is 

proportional to the energy loss in a square-network model for which M = 11 and /V = 21. This relationship, 
calculated from the data of table 3, is shown by the solid line. A relationship having the same value for 
8 = 0.10, but for which the ordinate is of the form K8 2 , where K is a constant, is shown by the dashed line. 
A corresponding relationship of the form K'8 1 ^ is given by the dotted line. This latter relationship provides a 
better fit to the data than does the quadratic relationship. 

The relationship between energy loss and local breaking force in an experimental test specimen is not 
known, and it does not seem feasible to determine it experimentally. It is possible that further studies with 
other more realistic mathematical models may provide more information on this relationship. 

5. Interpretation of the Energy Parameter. 

In the introduction an energy parameter E for characterizing the adhesion between pulp fibers was 
described. This energy parameter is obtained experimentally by elongating thin web-like specimens in a 
sensitive tensile tester, and measuring the energy loss incurred each time a bond breaks, as denoted by a 
sharp drop in the force-elongation curve. The energy losses for a selected number of bond breaks are then 
averaged. 
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The parameter E is supposed to provide a measure of the adhesive force between fibers, but the previous 
discussion shows that E depends on the structure of the specimen network as well. The meshes in a specimen 
network range over a variety of sizes, and the energy loss incurred by a break depends upon the size of the 
hole that opens up. If all the bonds in the neighborhood were of the same strength the parameter E would 
roughly measure that strength, but could still take on a range of values depending on the size distribution of 
the meshes in the specimen network. 

The force of adhesion between fibers in the network varies from bond to bond. It has just been shown that 
the energy loss incurred as the result of a bond break is not linearly related to the local breaking force at the 
bond. Therefore the average of a number of energy losses does not reflect a simple average of the local 
breaking forces. 

Results of calculations with square-network models suggest that the energy loss has only an approximate 
linear relationship to the mesh length. This introduces a slight additional complication in the interpretation of 
what the E parameter actually measures. 

The interpretation of the E parameter may perhaps be better expressed by the following formulation: Let the 
energy losses be measured for each of a series of n breaks, and let the energy loss AE, for the i ih break be 
given by 

AE,- = K t f t "g(xd 

Ki is a constant of the network that may vary from break to break if the network is significantly altered by the 
breaks./* is the local force on the bond at break, and a is a constant having a value probably between 1 and 
2. For a square-network model a has a value close to 1.5. g(xj) represents the functional dependence of AEj 
upon a characteristic mesh lengthy, associated with the bond i. This relationship is probably almost a linear 
one. 

Under these assumptions the characteristic energy E is given by 

1 ' = " 

E = -2 Ktfi a g(xi) (4.7) 

n /=1 

This formula gives an approximate idea of the nature of the parameter E. However one should remember that 
this formula is based on the results of square-network model calculations, and that the square-network model 
is only a crude representation of the paper networks that are tested experimentally. It is likely that work with 
more appropriate models will provide an improved understanding of the nature of E. 

6. An Alternative Force Drop Parameter. 

Adhesion between pulp fibers could be characterized alternatively by means of a force drop parameter F . 
This parameter could be obtained by averaging force drops incurred in a series of bond breaks when a test 
specimen is elongated. The parameter F has not been tested experimentally, but its use as a parameter to 
characterize adhesion seems feasible. 

In order to calculate an energy parameter E it is necessary to know the slope of the force-elongation curve 
in the vicinity of a force drop. The force drop parameter F could be calculated more simply, as only the force 
drops are used. 

In order that the force drop parameters be comparable, the test specimens must all have the same initial 
length. In other respects the force drop parameter F is similar to the energy parameter E. Both parameters are 
sensitive to mesh size and mesh size distribution in the test specimen network. Neither parameter is linearly 
related to the average of the local bond breaking force. However the functional dependence of F upon the 
average bond breaking force is different from that of E. 

According to eq (4.3) the local force drop resulting from a bond break is proportional to the quantity 
(1 — T)8. A plot of (1 — T)8 vs 8, which is proportional to the bond breaking force, is depicted in figure 8 for 
a square-network model with M = 11, N — 21. This curve can be fitted approximately by the curve 
(1 — T)8 ~ 0.00687 16 S 05 . It can be inferred from this result that the parameter F might not increase as 
rapidly with increase of local adhesion force as would be desired. Parameter E with its possible/ 1,5 
dependence may be superior in this respect. 
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The dependence of parameter F upon local adhesive forces and upon mesh size is somewhat like that of 
parameter E. A suggested formula to express this dependence would be similar to eq (4.7) with different 
values for K and a, and a slightly different quasilinear function g. 
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